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The effects of a large-scale shear on the energy spectrum of a passively advected scalar 
field are investigated. The shear is superimposed on a turbulent isotropic flow, yielding 
an Obukhov-Corrsin &~ 5 / 3 scalar spectrum at small scales. Shear effects appear at large 
scales, where a different, anisotropic behavior is observed. The scalar spectrum is shown 
to behave as /c _4//3 for a shear fixed in intensity and direction. For other types of shear 
characteristics, the slope is generally intermediate between the —5/3 Obukhov-Corrsin's 
and the —1 Batchelor's values. The physical mechanisms at the origin of this behaviour 
are illustrated in terms of the motion of Lagrangian particles. They provide an explana- 
tion to the scalar spectra shallow and dependent on the experimental conditions observed 
in shear flows at moderate Reynolds numbers. 



1. Introduction 

Two major regimes are known for a passive scalar field transported by a turbulent 
flow. The first, known as inertial-convective, refers to the range of scales where the direct 
effects of the fluid viscosity and the scalar diffusivity are both negligible. The predicted 
behaviour of the spectrum E e (k) of the scalar field 9 is | |Obukhov (1949]l|Corrsin (1951)| |: 



E e (k) = Ce- 1 ' 3 eg fc- 5/3 . (1.1) 

Here, k denotes the wavenumbers, C is an adimensional constant, e and eg are the average 
dissipation rates of the kinetic energy and the scalar variance, respectively. The second 
regime, dubbed viscous-convective, is at scales smaller than the former and holds for 
weakly diffusive (alternatively high Schmidt, or Prandtl number) scalars. The velocity 
field is now smoothed out by the molecular viscosity, linearizing the dependence of the 
velocity increments on the separation between the two points where the increment is 



taken. That yields the Batchclor behaviour (Batchelor (1959)) 

Eg{k) cx AT 1 . 



(1.2) 



Scalar spectra have been measured in a number of laboratory experiments. None of 
these pure power laws with the above anticipated exponents are routinely observed. The 
existence of a fc _1 in the viscous-convective subrange has been questioned based on high 
Reynolds and Schmidt numbers experiments ( Miller fc Dimotakis (1996)] Warhaft (2000)| 
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Yeung, Xu & Sreenivasan (2002)). As for larger scales, scalar spectra do behave compat- 
ibly with a power law, but their slope is often found to deviate from the Obukhov- 
Corrsin 5/3 value. Deviations are particularly important in turbulent shear flows, where 
the observations indicate that the spectra strongly depend on the Reynolds number Re 
(|Mestayer (1982)[|Sreenivasan (1991)]|Sreenivasan (1996)||Miller k Dimotakis (1996)[|Warhaft (2000) I 



and on the detailed set-up of the advecting turbulent flow. The prediction (|f .lfl is in par- 
ticular recovered for very large Reynolds numbers only (Re\ > 2000). Deviations are 
observed also when the velocity field displays a convincing fc -5 / 3 spectrum (see, e.g., 
breemvasan (1996)| |yillermaux, Innocenti fc Duplat (2001)| |. 



Our aim here is to investigate the physical origin of the behaviour experimentally 
observed in shear flows, briefly recalled in § 2. We specifically consider a model where the 
velocity is decomposed into a linear shear, superimposed onto a turbulent flow supposed 
to obey Kolmogorov's scaling. As in thermal convection (see, e.g., |Siggia (1994)] ), the 
presence of the shear introduces a new length-scale r c . The scalar spectrum crosses from 
the classical behavior (|1.1() at scales r <C r c to a different anisotropic behavior for r ^> r c . 
As discussed in § 3, Lagrangian dimensional arguments allow to predict the scaling of 
the scalar spectrum in the new shear-regime for a variety of shear properties. A —4/3 
scaling is, in particular, found for a shear fixed in intensity and direction. A test of 
the Lagrangian dimensional arguments is provided in § 4, where we consider a simpler 
turbulent flow with a short correlation time. The assumption is unphysical and the very 
same Lagrangian arguments of § 3 yield quantitatively different predictions; the point is 
however that these predictions can be stringently tested against the exact solution for the 
scalar spectrum which can be obtained in this special case. The final Section is devoted 
to discussions and conclusions. The Appendix is reserved to technical material used in 
§4. 



2. Scalar spectra in experiments 

Turbulent flows might be divided into two broad categories: shear flows (jets, wakes, 
and boundary layers) and flows where the effects of large-scale shears are tamed (e.g., 
grid-generated turbulence) . The latter are the laboratory flows providing the closest phys- 
ical realization of the homogeneous and isotropic turbulence which most of the theoretical 
fluid dynamics deals with. 

Various possible ways to inject scalar fields in turbulent flows have been employed. 
In grid turbulence, temperature fluctuations are generated by heating either the grid 
itself or small wires located thereafter. In shear flows, one possibility is to weakly heat 
the jet or the wall of the boundary layer. Alternatively, a colorant dye can be injected 
past the flow source, at a scale possibly different from that of the stirring source. Tem- 
perature time series are taken by a thermometer in a fixed location and, as customary 
in most fluid dynamics experiments, time distances are converted to space distances by 
Taylor hypothesis. Scalar spectra or, equivalently, second-order structure functions, are 
thus measured. The use of fluorescent dyes offers the possibility to extract cuts -usually 
two-dimensional- of the scalar field from the intensity field. Simultaneous velocity mea- 
surements are usually taken by anemometers. The reader interested in a detailed discus- 
sion of the experimental techniques is referred to Sreenivasan (1996)| |Warhaft (2000)1 
Villermaux, Innocenti fc Duplat (2001)| 



The observations gathered in the literature might be summarized as follows. In grid- 
generated turbulence, the prediction (|1.1() is well verified. That holds also at moderate 
Re\, when there is just a tiny range of scales where the spectrum of the velocity field 
displays a clean —5/3 slope. Some results obtained heating directly the grid used to 
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generate the turbulence have led to some criticisms (Warhaft (2000)) on account of pos- 
sible correlations between the temperature and the velocity. It seems, nevertheless, a 
well-established fact (see |Jayesh, Tong fc Warhaft (1994)[ |Mydlarski fc Warhaft (1998)| 



Warhaft (2000) I that the slope of the scalar spectrum increases with Re\ from values 
~ 1.5 to 5/3, reached at Reynolds numbers Re\ > 200. Small deviations observed at 
lower Re\ are confidently ascribed to finite Re -number effects. 

The situation in turbulent shear flows is more intricate. Experiments with different 
Re\'s but similar non-dimensional shear rates have been collected and compared by 



Sreenivasan (1996) The main observation is that the slopes of the passive scalar spectra 



display a strong dependence on the turbulent intensity. Their slope is approximately 1.3 
at low Re\ and attains the predicted 5/3 value at Re\ > 2000 only. Strong anisotropic 
effects are also observed for the velocity energy spectrum. Indeed, its longitudinal com- 
ponent reaches the expected —5/3 behaviour at moderate Re\(~ 50). Conversely, its 
transverse component behaves more similarly to the scalar field, with a slope continu- 
ously increasing with Rex and saturating to 5/3 only for Re\ > 3000. The experiments 
by Villermaux, Innocenti & Duplat (2001) report even more severe deviations from p.ljl . 
with a slope close to the Batchelor unit value, in spite of the —5/3 spectrum observed 
for the velocity field. The scalar injection scale was, although still lying in the inertial 
range of scales, smaller than the turbulence forcing scale in that case. 



3. Shear effects on Lagrangian dynamics 

The experimental observations presented in the previous Section naturally lead to 
surmise that the deviations observed in the scalar spectra have, among other possibles 
causes, their origin in the presence of a sustained large-scale shear. The interplay between 



diffusion and shear is a classical problem first examined by Leveque (1928) for heat 
transport across boundary layers. Qualitatively, the flow stretches the distances between 
particles along the shear and sharpens the scalar gradient perpendicular to it, therefore 
altering the usual diffusion law. 

In order to quantitatively assess the impact of shear on scalar spectra, we introduce 
the following simple model for the velocity field v(r,t): 

v(r,t) — o-yx + u(r,t). (3-1) 

The flow is thus the superposition of an average linear shear of intensity a directed along 
the x— direction and a turbulent fluctuating field it, supposed to obey Kolmogorov's 
scaling, for the sake of simplicity. The spectrum of the velocity field obeys in particular 
the classical fc _5//3 law. 

A passive scalar field, 6(r, t), transported by the above flow is governed by the standard 
advection-diffusion equation: 

d t 6 + ayd x 6 + u ■ V9 = kA9 + f , (3.2) 

where k is the molecular diffusivity and / is an external source of scalar fluctuations. 
The source is needed to maintain the system in a statistically steady state and we shall 
suppose that the forcing has a characteristic length Lf. 

The transport equation l|3.2|l is equivalently recast in terms of the Lagrangian trajec- 
tories p(t) of tracer particles. They obey the stochastic differential equation 

dp = v(p,t)dt + V2KdW, (3.3) 

where W is an isotropic Brownian motion and v is given by 13.1(1 . The equation for 
the scalar field along the trajectories is the ordinary differential equation dO/dt = /, 
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easily integrated. The scalar second-order correlation function of the scalar field is then 



expressed as (see, e.g., Falkovich, Gawedzki & Vergassola (2001)1 



C 2 (r,t) = (6(r,t)e(0,t)) = [ ds f dRp(R,s\r,t) X (R) ■ (3.4) 

Here, R is the separation between two tracer particles and p{R, s|r, t) is the propagator, 
i.e. the probability, averaged over the realizations of the velocity ensemble and of the 
Brownian noise, that a pair of particles is found separated by R at time s, conditional 
to its separation r at time t. The function x(i?) having its support at scales R < Lf, 
the integral (|3.4[1 is roughly equal to x(0)7 1 (r, L/), where T(r,Lf) is the residence time 
at distances < Lf for a pair of particles initially separated by a distance r. Analogously, 
the second-order structure function S 2 (r) = 2 [C 2 (0) - C 2 (V)] oc T(0,L f ) -T(r,L f ). 
The incompressibility of the flow ensures that the scaling of the latter quantity can be 
estimated as the time for two particles, initially coinciding, to reach a separation r. 

The problem of determining the scaling behaviour of scalar spectra is thus reduced to 
the study of the evolution of the separation between a pair of particles. 

For the case without shear, a — 0, classical dimensional arguments give the Richardson 
law: (R 2 (t)) ~ t 3 , which translates into a 2/3 scaling exponent for S 2 (r) and the ensuing 
Obukhov-Corrsin —5/3 slope for the scalar spectrum. 

Let us now consider the case with shear, a ^ 0. The components of the separation 
between two fluid particles obey the equations of motion: 

R x =<jR y +6u x (R,t), (3.5) 
R ± =Su ± (R,t), (3.6) 

where the dot indicates the time-derivative, 5u is the velocity difference between the two 
particles and R± indicates the components of the velocity orthogonal to the shear. At 
small enough times, the turbulent component, which scales with exponent 1/3, dominates 
over the linear shear and Richardson's behaviour holds. The crossover to a different 
anisotropic behaviour occurs at r c ~ er -3 / 2 . Those are the separations where the shear 
and the turbulent components become comparable. The time required for two particles 
initially coinciding to reach those scales behaves as t c oc rr -1 . For t ^> t c , the form of the 
Lagrangian equations naturally suggests that the parallel and the transverse components 
scale differently with time: 

(R 2 X ) m (R 2 ) ex t 2a , and {R 2 ± ) ~ t 2/3 . (3.7) 

Inserting this ansatz into (|3.5|l and l|3.6|l , neglecting the Su x term in l|3.5() and using the 
Kolmogorov's scaling 5u±(R,t) oc i? 1 / 3 ~ W 3 , we obtain the relations: a — 1 = and 
(3 — 1 = a/3. We thus end up with the predictions a — 3, (3 — 2 for the Lagrangian 
separations, that give 

CsWocr 1 / 3 ; E(k) oc fc~ 4 / 3 , (3.8) 

for the correlations of the scalar field, by using the arguments following (|3.4|l . 

It is worth emphasizing that the scalar statistics is sensitive to the shear character- 
istics, as it should be expected for any passive transport. For example, if the direction 
of the shear were to rotate rapidly, isotropy would be recovered and the stretching of 
the separations among the particles would be dominated by the linear shear component, 
leading to an exponential-in-time separation and a fc" 1 behaviour. The physics is the 
same as in the Batchelor regime, but the behaviour now holds at large scales. The rapid- 
ity of the rotation should be gauged with respect to the Lagrangian turn-over time: if 

3/2 

denotes the typical time of rotation, the fixed limit (|3.1|l holds for scales r <§; rJ and 
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the rapidly rotating asymptotics in the opposite limit. As for the fluctuations of the shear 
intensity, let us similarly denote by T a its correlation time. The intensity is effectively 
fixed in time, as in (|3.1[1 . for scales r <C t^ 2 . In the opposite limit, the intensity might 
be taken as a random process with a short correlation time (a(t)a(t')) = a 2 8{t — t'). 
Inserting the ansatz (|3.7|) into 1)3. 5JI leads now to the relation a — 1 = — 1/2 + /?, giving 
a slope —13/9 for the scalar spectrum. The geometry of the shear and the ratio between 
its time-scale and the turbulent turn-over times are thus seen to affect the properties of 
a passively transported field. The point of interest to interpret the experimental data is 
that the exponents of the spectra are in all cases smaller than 5/3, implying that the 
presence of a shear generally induces a crossover to a shallower spectrum at the large 
scales. 



4. A solvable case 

The aim of this Section is to test the dimensional arguments presented in the previous 
section by investigating the simpler case of a turbulent flow u belonging to the Kraichnan 
ensemble (Kraichnan (1968)1. The rationale is as follows. The short correlation time of 



the velocity is unphysical and the Lagrangian predictions (|3.8fl are modified thereby (the 
scaling —4/3 of the scalar spectrum becoming, for example, —11/9). Yet, the point is 
that one can compare those predictions versus the solution of the exact equation for the 
scalar spectrum. The advantage over a direct numerical simulation of the more realistic 
flows in Section 3 is that the range of scales available here is much more extended and 
the comparison is therefore quantitatively more stringent, albeit less direct. 

The field u in 1)3. l|l is taken now as an incompressible, statistically isotropic and ho- 
mogeneous Gaussian field of zero mean and correlation function: 

([« i (r,t)-« i (0,t)][u i (r,0)-u i (0,0)])=*(t){Dr« [(1+0%-^]} • (4.1) 

Here, we specialize to the two-dimensional case for the sake of simplicity: this choice bears 
no qualitative consequences on the results. The parameter D measures the turbulent 
intensity and the velocity field is assumed to be scale-invariant with roughness exponent 
£ € [0,2]. The presence of an ultraviolet viscous cutoff mimicking the Kolmogorov scale 
might be considered, but it will not be needed here. 



Using standard methods of Gaussian integration by parts (see Falkovich, Gawedzki & Vergassola (2001) I 
we can re-adapt the Lagrangian arguments of the previous section to the case of a Kraich- 
nan flow. The Lagrangian separation law for the isotropic case is: (R 2 (t)) oc t 2 /( 2 ~?) . Note 
that the Richardson scaling is obtained for £ = 4/3, differing from the 2/3 of real tur- 
bulent flows due to the (5-correlation in time of the velocity field. For the anisotropic 
case of a shear fixed in intensity and direction, the relations between the exponents a 
and /? stemming from (|3.5I) and l|3.6l) read now: a — 1 = j3 and 2/3 — 1 = 2a£. Solving 
these relations gives the following behaviors of the Lagrangian separations and the scalar 
spectra: 

(R 2 X ) « (R 2 ) oct 6 /^); (i^ t 2(i+0/(2-S). K fc -(B-e)/3 _ (4 2) 

For £ = 4/3, (|4.2|l gives a —11/9 exponent for the scaling of the scalar spectrum. 

The advantage of Kraichnan's flows is that the behavior of the second-order correlation 
function Ci(r) at the stationary state can be determined exactly and used to test the 
previous dimensional arguments. In the shear-free case, a = 0, the problem is isotropic 
and the expression of C2(r) is well-known analytically. The inertial-range scaling be- 
haviour for the second-order structure function S^M = (2- < '°)_D r2 ~^' where x(0) is the 
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Figure 1. Scalar spectra in the x— direction, Eg(k, 0), obtained by fixing D = 0.05 and varying 
a (from top to bottom, a — , 0.0004 , 0.0036 , 0.025). In the numerical computations Lt ox and 
Lf have been fixed to 128 3 and 10 5 , respectively - see the Appendix for the definition of the 
parameters - and the scaling exponent of the velocity £ = 4/3. The inset on the top right shows 
the rescaled spectra Eg(kr c ,0) vs kr c . At large wavenumbers, the fc _5//3 spectrum is recovered, 
while at small k the power law fc -11 / 9 , predicted by Lagrangian dimensional arguments, appears. 
The inset on the bottom left shows the compensated spectra Eg(k,Q)k 11 ^ 9 obtained by fixing 
a = 0.025 and varying L f = 10 4 , 10 5 , 10 6 , 10 7 . The number of harmonics t, the box size Lf, ox 
and the number of grid points have been varied accordingly. Note that the transition region to 
the asymptotic slope —11/9 is quite broad. 



scalar variance injection rate. The exponent 2 — £ is in agreement with the Lagrangian 
dimensional estimate. 

In the presence of shear, the equation for the second-order correlation function is still 
closed, but the problem is not isotropic anymore. Following the same procedure as in the 



shear-free case (Kraichnan (1968) I, one obtains 



dtC-i + a(r sin <j> cos <f> d r — sm 2 (j)d^)C2 = -d r (Dr 1+ ^ + 2n,r)d r C2 



r 

+ i[(l+£)Drf + 2K]^C 2 + X (r), (4.3) 

where r = (r, cf>) denotes the relative distance between the two points. The forcing is 
taken Gaussian, statistically homogeneous and isotropic, with zero mean and correlation 
function (f(r, i)f(0, 0)) = 5(t)x(r). In the simulations, the correlation function of the 
forcing is taken as x( r ) = exp(— r 2 /2Lj). Contrary to the <5-correlation of the velocity, the 
previous hypotheses on the forcing are not restrictive and might be easily relaxed. The 
anisotropy induced by the shear makes that C2 will depend both on r and on the angle 
<fi between r and the x— axis. The steady-state solution of Eq. Ij4.3|l cannot be obtained 
analytically anymore and we have to resort to numerical methods (sec the Appendix for 
details). 

The existence of a crossover scale, r c , separating the scales dominated by the turbulence 

and those affected by the shear, is easily recognized by a direct inspection of Eq. I|4.3|l . A 

1 

simple balance of the shear and the eddy-diffusivity terms in Eq. 1)4.3(1 gives: r c ~ (— ) 2-4 . 
The crossover is well evident in Fig. 1, where the scalar spectrum for £ = 4/3 in a wide 




Figure 2. (a) The scalar second-order structure function measured in the a;— direction, ^(r, 0), 
fixing a and varying Lj (as in the bottom inset of Fig. 1). The inset shows the local slopes. 
Note that the approach to the asymptotic scaling r ' 22 is even slower than in the spectra, (b) 
The scalar second-order structure function measured at different angles, S2(r,<f>), for the case 
a — 0.025 and Lf = 10 7 . The statistics becomes more and more anisotropic as larger scales are 
considered. A similar behaviour is observed for the spectra (not shown). 



span of scales is reported. The results are in agreement with the Lagrangian predictions 
l|4.2jl . with E e {k) ~ fc~ 5/3 , for kr c > 1, and E e (k) ~ k~ n / 9 asymptotically reached at 
the large scales. Moreover, as shown in the top right inset, the spectra obtained with 
different shear rates all collapse by the rescaling of the wavenumbers kr c . 

The corresponding behaviour in physical space for the second-order structure function 
is shown in Fig. 2a. Two regimes are again clearly visible. A point to be remarked in 
the inset is that the convergence to the asymptotic slope, 62(7*) ~ r 2 / 9 , is even slower 
than in fc-space. The behaviour of the second-order structure function for different ori- 
entations is shown in Fig. 2b. As expected in a range of scales affected by the shear, the 
anisotropy is strong and the crossover between the two previous regimes depends on the 
angle of measurement. It would be very informative to have similar measurements of the 
anisotropy in experimental conditions. 

Note, as a side remark, that an asymptotic expansion valid for small angles can be 
found analytically. Indeed, the correlation function C2 can be expressed as 

C 2 (r x: r y ) = (D/ C r)-V3 r (2-e)/3 /[rj//(i)r i + e /cr) i/3 ] ? (44) 
where the angular part / obeys the equation: 

For small angles, i.e. w < 1, the operator appearing in 1(4. 3|) reduces to <Jr y d Tlc — D(l + 
£,) r idr y an d it is easy to check that the solution to 14.5J1 has the form 14.4J1 with the 

function / behaving as f(w) = -1 - ^TTo ~ eTTTf) +°( wi )- 

An even closer contact to the Lagrangian arguments 14.2(1 is made by directly simulat- 
ing the equations of motion for the tracer particles. The method used for the simulations 



is the same as in Frisch et al. (1998) Fig. 3 presents the evolution of the particle sep- 

aration for £ = 4/3. As long as (R 2 (t)) < r 2 ~ (— j 3 ^, the Richardson prediction is 
recovered and the anisotropy induced by the shear is negligible. As the separation be- 
comes larger and larger compared to r c , the trajectories are more and more affected by 
the shear and (R 2 ) « (i? 2 ) » In this range of scales, (i? 2 ) and (R 2 ) grow with 

two different power laws and the local slopes in the inset of Fig. 3 are found to be in 
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Figure 3. The curves of the Lagrangian separations as a function of time, for a Kraichnan 
flow with exponent £ = 4/3. At small times, {R 2 x {t)) = {R 2 y {t)) = (R 2 (t))/2 ~ t 3 , as expected 
according to the Richardson law. At large times, (R 2 (t)} ~ {R 2 (t)} ~ t 9 and (Ry(t)} ~ t 7 , in 
agreement with H4.2I . The inset displays the local slopes as a function of time. Note the wide 
crossover region between the two regimes. The numerical integrations of Eqs . 13.51 and 13.61 
have been performed following the point-splitting scheme described in |Frisch et al. (1998)| i.e. 
with molecular diffusivity k — and taking the initial separation small but finite. 



agreement with (14.2(1 . The predictions ((4.2(1 and the numerical simulations are found to 
be in agreement for values of £ ^4/3 as well (data not shown). 

An interesting remark is that the Lagrangian numerics evidences a wide crossover 
region between the two different regimes (see Fig. 3). That is the reason for the very slow 
convergence of the scalar structure functions and spectra to their asymptotic slopes. In 
general, detecting the transition between the inertial-convective and the shear regimes 
demands a very broad range of accessible scales. 

We finally consider the case where the shear is random. If its direction is fixed and 
the time-scale of its intensity is short, the resulting exponent of the scalar spectrum is 
— 2+£/2. If both the direction and the intensity of the shear rapidly fluctuate in time, i.e. 
the shear is a Kraichnan flow with £ = 2, isotropy is recovered and the scalar correlation 
function in the stationary state obeys (l/r)<9 r (or 3 + Dr 1+ ^ + 2nr)d r C2 = —x( r )- The 
equation is easily solved analytically and the ensuing slopes of the scalar spectra are 
— 3 + £ at small scales (large wavenumbers) and —1 at large scales (small wavenumbers). 



5. Conclusions 

The results presented here illustrate the presence of a new shear regime for a passive 
scalar transported in turbulent shear flows. The scale of crossover between the classical 
Obukhov-Corrsin range and the larger scales, where the effects of the shear are felt, is pro- 
portional to the ratio between the turbulent and the shear intensities. The range of scales 
where the shear-regime is observed becomes then parametrically large with increasing the 
shear rate. In experiments where the latter is fixed, increasing the turbulence intensity 
has the double effect of decreasing the Kolmogorov scale and increasing the crossover 
length r c . The Obukhov-Corrsin scaling is therefore observed over a broader and broader 
range of scales. When r c is much larger than the integral scale L of the scalar, the picture 
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at scales < L becomes essentially shear-free. For r c ~ L, the presence of the shear-range 
is still relevant due to the crossover effects which induce a scalar spectrum shallower than 
—5/3 and anisotropic. Those effects disappear quite slowly with the Reynolds number 
due to the broad span of scales required for the crossover. Furthermore, we have shown 
that the scalar spectrum in the shear-region is sensitive to the geometry and the fluc- 
tuations of the shear. Experiments with different set-up's, where the shear components 
have drastically differing properties, might then yield different scalings at asymptotically 
large scales. All are shallower than —5/3, though, and therefore inducing crossover effects 
that tend to flatten out the scalar spectrum. In those conditions where shear fluctuations 
in intensity and direction are kept controlled, we specifically predict the existence of a 
fc -4 / 3 asymptotic regime. The previous picture seems to account for the major points in 
the experimental observations presented in § 2. 

Two issues for future investigation are the following. Theoretically, it is quite likely that 
the shear will affect the scalar higher-order statistics; the resulting behaviours remain 
to be clarified. Experimentally, it would be of major interest to have a set-up where 
the Reynolds number is kept fixed and the shear rate is increased, in well-controlled 
conditions, so as to enlarge the shear-regime and quantitatively analyze its properties. 

The authors acknowledge the hospitality of the Erwin Schrodinger Institute in Vi- 
enna. This work has been partially supported by the EU under the contract "Non-Ideal 
Turbulence" HPRN-CT-2000-00I62 and by Cofin 2001 (prot. 2001023848). 



Appendix A. Equations for the correlation function and numerical 
computation 

The anisotropy introduced by the shear can be handled by using polar coordinates, 
which allow to decompose the correlation function as: 

oc oc 

C 2 {r, t) = C 2 (r, cb, t) = F Q {r, t) + F t (r, t) cos(2£</>) + ^ G e (r, t) sin(2ty) . (A 1) 

i=i e=i 

Here, i labels the order of the harmonics. Note that only even harmonics are present 
for symmetry reasons. Substituting IjA 1(1 into 14.3(1 . one obtains the following set of 1+1 
dimensional partial differential equations: 

dfFi +a r -[d r G t+1 - drGt-i] + l)Gt+i - 2tG t + {£- l)G<_i] = 

1 M 2 

-drpKr + Dr^drFe-pK+il + ODr^—Fe + Sevx, (A 2) 

d t G e -a r -[d r F l+l - (1 + S^drFt-!] --[(*+ - + (£- - 

1 AC 2 

-d r [2kt + L»r 1+ «] d r G e - [2k + (1 + £)i>«] — G t . (A3) 

Note that the coupling among the harmonics is indeed due to the shear term. 

Those equations have been numerically integrated for times long enough to reach a 
stationary state. The correlation function Ci(r, <fi) is then computed at different ori- 
entations (j> from Fi{r) and Gt(r) by using ((A 1() . The second-order structure function 
S2{r,4 l ) = 2(C2(0) — Ci(r, (j))) easily follows, as well as the scalar spectrum Eg(k,4>), 
which is obtained by a one-dimensional Fourier transform. 

The equations have been solved for various choices of the parameters. In practice, we 
have fixed D and varied a and Lf. As for the boundary conditions at the origin, we have 
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taken F£(0) = 0, G e (0) = for each I and F e (0) = for I > 0. The infrared boundary 
condition has been fixed by setting to zero the values of all the functions Fg{r) and Ge(r) 
beyond r = L\, ox . The maximum size of the system, Lb ox , has been taken much larger 
than the forcing scale Lf to properly resolve the anisotropic contributions. 
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